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ABSTRACT 


"ft 


/ 


A new class of problems has been analyzed to 

f?ed^inite U sefSf°ho^enloS rotlrilis 

and discussed using specific examples. 

INTRODUCTION 

The problem of optimizing strength properties of inhomogene- 
ous targets under impact of tapere „,]???* lve criteria of optimal 
considered in the developed in thi^and^h^f oil owing studies Cfor 

anal yt i cal rel at ion between material hardness and density a 

°^ SS In'" the present°investigation°the range of control (unctions 

nf^SidleMmtioKs ?41 when Estimating the 9 ®necessary optimiza- 

?10n analysis 1 of 9 pene?ration"is Sffi^th. model 

of radial nole extension [51. 

ANALYSIS 

1. Penetration model. 

indenter. rPn + k 0 nrp<?«;iire actino on the indenter within 

the dtsUnoe n V= x-?' surface can be written 

as 


p = a/ 2 Ppu 2 {far/a?; 2 [ 8 f £ >£/a+^] + ra 2 r/a? 2 } + 

+ (l/25pu QCcjrdr/dt, + f !/25d s [i + SCc^ 


Cl. 15 
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there r = rC?) is the expression for the generating line of the ^ 
axisymmetric indenter, p is the density, SCe) -Incite), e = 
L/icr^l+vJi , v is the indenter current velocity, cr is the yield 

stress, v is Poisson ratio, £ is Young's modulus. 

form " ie e£ I ua *'* on m °tion for the indenter of mass M has the 

C l/2)Md(v z )/dx = -2n f pCprC^Xdr/dV^ Cl. 2) 

with the initial condition uCO) = v° 

0 

2. Optimization problem. 

It. is convenient to represent the penetration process as a 
system of differential equations relative to a vector of Dhase 
coordinates y with _u as a control function H 

dy/dx = J Cy,iD, yCO) = y o (2.1) 

and to define Freshet’s differentiable functionals by 

F tvC.O y bl = f pCxOdx, (2.2) 


(2. 1) 


- J pCxjdx, (2.2) 

0 

F j [u(.l y b] = y 1 (b) = 0 (2.3) 

, -iH^ertion of additional phase coordinates and use of (1 1) 
and (1.2) reduces (2.1) to 

dyVdx = ~2(y l y 2 + df)/(d + y 4 J>, dy 2 /dx = y s - ACx-hO-h, 
dy 3 /dx =y e -B(x-hO'h y dyVdx =2y 7 -CCx-h>h 2 , dyVdx =A(x)-A(x-hO , 
dy e /dx =B(xO-B(x-hO t dyVdx =y e -C(x-hO-h y dy B /dx =C(x)-C(x-h ) , 

y l CO) = v z Q , y *(00=0, Ci = 2 y 3 y . . . 8) , (2.4) 

where d = Ctg a) -2 , d g = M dVrc, y 1 = u 2 , 

ACx) - 0 Cif x<0), pCx) [bCcCx)) - e(x)/(l+£Cx )) j (if x>0 J, 

BCxl = 0 (if x<0) y <? s (xO U + B(e(x^ ] (if x>01 y 
C(x) = 0 (if x<CO, pCxJ 9(e(x)) (if x>0) y 

The size, number and class of materials to be used in the 
tiel 6 t layers are s P ecified by a distribution of material proper- 

uC *> = <u s : xcfx s> x s+i ;, s=T7n > y x = O y x n+i = b y 
where n is the number of layers. The value of u £ "belongs to a 
finite set U which corresponds to a given set of materials u Q U = 
(U x ,U z . . . to. Here u & is the material in the s th layer, U is the 
material number and q is the material quantity. 
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r The stated optimization problem suggests that from all piece-1 
wise continuous functions u(x) 9U and numbers b >0 one should 
choose a control iu° Cxi, b°> which will provide minimum for the 
functional C2.2) under the limiting conditions (2.1), (2.3). The 
quality criterion may be refered to as a specific plate mass (2.2) 
subject to u(b) = 0 (under the requirement of arrested indenter). 


3. Necessary conditions of optimization. 

A discrete character of the control function range doesn’t 
allow to generate small variations in the norm lldull =max |<5u|. The 

xS [ o, bJ 

disturbed control may be written in the form 


ut x) = 


- f W . 

- w°. 


X971, w£U 

x$m.> 


(3.1) 


where m£ [0,b°] is the set of measure zero. 

An equation for the system (2.1) is expressed in terms of 
variations and the main terms of functional increments (2.2), 

(2.3) are given by „ 

- <5Cdy/dx) - df/dy 6y = f(y,iD - f(y,u° 


6F - f [pCw3 - pCu°)]dx + p[u°Cb° 3J6b ° , 


(3.2) 


m 


6F i = Sy'Cb 0 ) + 7[y(b°),u°](55 c 


Using the Lagrange identity and desired limiting conditions 
for the disturbed trajectory one finds an expression for <5fc> 


6b c 


= [f// t fy(b°) ,u°]J J y/f/Ey ,iD - f[y,u 0 )ldx (3.3) 


(3.4) 

(3.5) 


where the conjugate vector-function. yj_ satisfies 

dy/dx = - d]/dy y 

Variation of the minimized functional is written as 

6F Q iu(.hb] = J [H(y } y,u°) - H(y,y,w)]dx 
m. 

In order to make the control function optimal it is necessary 
to follow the principle of maximum 

H(y t y,u°) = max HC y,y/,iiD (3.6) 

w^U 

An expression for M is given as 
H - D^p s Cx)ACw3 + y e Cx)BCw) + y^xlCCwiypCw), m£[b°-h t b° J, 

H - f^ s (x)-y z (x+hJh-y^(x+hj'^A(w)+ ^y^Cxl-y^x+hJh-y^Cx+h^BCw). . . 


+ ^y^xJi-y^x+hJ-y^Cx+bJh-y^Cx+tOh 2 j j- - pCw m£tO,b°-h] , (3.7) 
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where D = p[u°Cb°)]/f f yC 5°),u°3, xQn. 1 

4. Geometrical interpretation and qualitative conclusions [61 . 

n 

CD Function M can be expressed as H = ^ p . Cy , yD<p i C w ) , <p i are 

continuous functions of x. The function H given in p is refered to 

as a hyperplane of support to a vector-gradient, which defines 
direction of increase grad H - <ju , . . . ,ju n >. From this follows 

that M approaches maximum at one of the vertices of convex poly- 
hedron 0. , which represents a convex shell of the point set 
PjCw ) , s = 1 ,q . The remaining materials of the set U can 

be excluded from a further discussion. 

Cii) The continuity of p Cx) implies that at any vertex of 

the polyhedron Q there is a hypercone K the interior of which may 
contain qrad H at slight variation in x? Cx ,x )c to, 53 and allow 
to satisfy the maximum condition. Thus, the optimal plate structu- 
re includes the finite number of layers of finite thickness. 

C i i it) Substitution of materials is expected to take place at 
the contact points x of the hypercone K and one of the polyhedron 
edges. It is to be noted here that the immediately adjacent mate- 
rials may be only there which match the adjoininq vertices of the 
polyhedron. 

Civ) It can be shown that from the entire set of materials 
assumed in the vicinity of the rear surface the preference should 
be given to material with minimal density. 

5. Numerical algorithm. 

Numerical procedure requires insertion of some admissible 
control function uCx) Q U and a small parameter x which describes 
the set of measure zero. Computational algorithm involves the 
uniform mesh x having the mesh spacing %■ The values of y and y 

are calculated at points x s + */<? and assumed constant for the 

segment ^ x s - x s + 1 ^- 

Solution includes the following steps: 

Ci) The system C2.4) is integrated and yCx) and b are defined 
at mesh nodes. 

Cii) Boundary conditions for conjugate functions are prescri- 
bed at the point x = b and the system C3.4) is solved. 

Ciii) A new value of u* on the segment m, is derived from 

condition HC . ,. ,u*) = max HC . ,. ,iD; if u* = u(x+ */<?.) this step 

s w^U s s 

is repeated for s = s+1 ; otherwise, a new control function is 
assumed u - u C if xQn ) , uCx 0 Cif x$m. 3 and calculation returns 

s J s s 

L J 
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f to step C i) . 

The procedure of improving control function proceeds like 
this up to the terminal point on the right of the interval. The 

process is completed as soon as uCxO remains constant at any s. 

RESULTS 

A set of materials contain annealed aluminum (a), aluminum 
alloy B-95 (b) , annealed titanium (c), titanium alloy BT-6 
(d), steel Ct. 6 (e), steel 12X2H4A (f). Material properties are 
given in Table 1. 

Table 1. 


Material 

Density 
p, g/sm 3 

Young’s 

modulus E, GPa 

Yield 

stress cr , GPa 

S 

Poisson 
ratio, v 

a 

2.8 

70 

0.06 

0.33 

b 

2.8 

70 

0.45 

0. 32 

c 

4.5 

110 

0.08 

0. 30 

d 

4.5 

120 

0.83 

0.32 

e 

7.8 

200 

0.21 

0.28 

f 

7.8 

200 

0. 11 

0.30 


For the case of dynamic penetration of tapered cylinder the 
optimal plate will consist of two materials (d)+(b). The relative 
front layer thickness of the optimal plate increases with the 
increase in a half-angle of the cylinder opening a and initial 
impact velocity v . 

Fig. 1 shows the decreased mass optimally T q of homogeneous 

plates made up of CiO, Cd) and (/) -materials as compared to the 
optimal one. Disadvantage of heavy materials (d) and (/) decreases 
with the increase of v q , since their fraction in the optimal plate 

is growing high. At v o <600 m/s the preference is given to a homo- 
geneous material (b) rather then C d) , while at v >600 m/s the 

O 

prefered material is of d- type. 

The results of present investigation agree qualitatively and 
quantitatively with data reported in C21. If instead of material 
Cf) one uses steel with the yield stress 1.5 GPa an optimal plate 
will consist of three layers. In case of a large choice of materi- 
als an optimal target structure will be multilayer. However, the 
main qualitative characteristic - a decrease of density and hard- 
ness with a distance from the upper to lower surfaces of the 
target - remains uncharged. 

It is to be noted here, that the usefulness of a soft rear 
layer in a target has been already justified but only in context 
of fracture behavior of material. From mechanical point of view 
the optimality of target structure predicted in [2] and in present 
investigation implies its high resistance to penetration while 
preserving the same ductile type of cratering. 
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